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Abstract 

We discuss the shape dependence of the finite-size scaling hmit in a strongly ani- 
sotropic 0{N) model in the large-A^ limit. We show that scaling is observed even if 
an incorrect value for the anisotropy exponent is considered. However, the related 
exponents may only be effective ones, differing from the correct critical exponents of 
the model. We discuss the implications of our results for numerical finite-size scaling 
studies of strongly anisotropic systems. 
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1 Introduction 



Finite-size scaling (FSS) ^EllSlll] is a very powerful tool that allows to extract infor- 
mation on the critical behavior of a system — which in principle can only be observed 
in the infinite- volume limit — from finite- volume results. In particular, the most recent 
Monte Carlo studies heavily rely on FSS for the determination of critical properties 
(see, e.g., Refs. [5ll^ll7|l§l l^inpiipi2pi3j for recent applications to the A^- vector model in 
two and three dimensions; the list is of course far from being exhaustive). For isotropic 
systems the recipe is well known. For instance, in two dimensions one considers a 
domain L x M and the finite-size scaling limit (,l.m oo, L, M ^ oo at fixed ^l,m/L 
and aspect ratio L/M. Here ^l,m is a suitably defined finite- volume correlation length. 
In this limit, long-distance quantities, e.g. the susceptibility, show a scaling behavior. 
For instance, if O diverges in the thermodynamic limit as Itl"^*^ for t = T — Tc — > 
(T is the temperature and Tc its critical value), then in the FSS limit one finds 

L, M) ^ L^'^/-fo{iLMt)/L, L/M), (1) 

where foi^iU) is a universal function. 

One may also ask what happens if one considers a different limit: ^l,m — *■ oo, 
L,M ^ oo keeping fixed the ratios ^lm/L and L/M^+^ with 5/0. If (5 < 0, M 
increases faster than L and it is easy to guess that we will obtain an effective strip 
geometry so that 

0{t, L, M) ^ L-^/''fo{^L,Mit)/L, 0). (2) 

On the other hand, if (5 > 0, M increases slowly, and, if we are able to keep ^,l,m/L 
fixed, i.e. to use L as reference box size (it is not obvious that this is possible), the 
domain effectively shrinks and becomes one-dimensional. Then, the question is whether 
a scaling behavior is still observed. One may imagine that the scaling function to be 
used is that of a one-dimensional system, but in this case it is unclear which exponent 
should be used in the prefactor. 

This problem may appear academic at first but it is motivated by systems that are 
strongly anisotropic. For these systems one introduces an anisotropy exponent A and 
the canonical FSS limit is obtained by taking ^l,m oo, L,M oo keeping fixed 
Cl,m/L and the aspect ratio L/M^~^^. Examples of such systems are provided by driven 
systems that admit nonequilibrium stationary states with strong anisotropy jl41ll5j . 
surface-growth processes ^Hl, Lifshitz points ^3^1^11213) uniaxial magnets with 
dipolar interactions j21ll22j . just to mention a few of them. A general approach to 
scale invariance in infinite volume for these anisotropic systems has been developed 
by Henkel [^. However, FSS is still poorly understood. An exact computation 
on a dimer model which undergoes an anisotropic phase transition was performed by 
Bhattacharjee and Nagle [21], while, in the context of the FSS for driven diffusive 
systems, a phenomenological FSS theory was proposed by Binder and Wang [2^] and 
by Leung I2(ip27j who also gave some heuristic arguments on the consequences of taking 
the FSS limit with an incorrect aspect ratio. 

Often in experimental or Monte Carlo applications the exponent A is not known 
a priori, and thus two questions naturally arise. First, if we consider the FSS limit 
with fixed L/M^+^+\ 5 / 0, do we still observe scaling? and if yes, with which 
exponents? Second, how do we determine A? An answer to these two questions is of 
utmost practical importance. 
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In this work we analyze the finite-size scahng behavior of the A'^- vector model in the 
N ^ oo limit, since in this case it is exactly solvable, see, e.g., Ref. |2SI- This model has 
been a classical example in FSS investigations, starting from the work of Brezin j29j . 
recently extended to several generalizations ^| , including long-range interactions and 
boundary effects. Here, we consider a general class of 0(oo) models which includes the 
classical isotropic short-range and long-range cases but also models in which the spin- 
spin coupling decays with different power-laws in different lattice directions, giving rise 
to a strongly anisotropic phase transition. For these models we analyze the FSS limit 
using finite boxes with arbitrary shape. 

We show that A is uniquely determined if one properly measures the correlation 
length. On the other hand, if only zero-momentum quantities are available, for instance 
the susceptibility X) A cannot be easily determined. Indeed, even if the aspect ratio 
is incorrect, i.e. 6^0, one still observes scaling. In some cases, it is even possible 
to observe two different FSS limits with the same data, one corresponding to the 
layer /strip geometry, the second one corresponding to a lower-dimensional system. 

Our analysis will be limited to systems below the upper critical dimension and 
above the lower critical one. Also, we will not address the question of the FSS limit 
at fixed vanishing magnetization, which is of relevance for lattice-gas studies and has 
already been discussed for isotropic systems in Ref. ^Tj . 

The paper is organized as follows. In Sec. [21 we define the model, the basic observ- 
ables, and discuss the large- limit. In Sec. 01 we derive the anisotropy exponent A 
that defines the aspect ratio and report the FSS functions for the susceptibility and 
the correlation lengths. In Sec. jll we discuss the noncanonical FSS limit in which the 
ratio M/L^+^+\ (5 / 0, is kept fixed as the size of the lattice is increased. The results 
of these two sections are derived in the Appendices. In Sec. [51 we present a simple 
numerical example: we consider the standard isotropic model with short-range inter- 
actions on a cubic lattice M x and show that two different scaling behaviors can be 
observed by keeping fixed the ratio M/L^/'^, in agreement with the theoretical results 
of Sec. 01 Finally, in Sec. [SI we present our conclusions and discuss the implications for 
numerical studies. In particular, we discuss how one can determine numerically the 
anisotropy exponent A. 



2 The model 

We consider a (i-dimensional hypercubic lattice Z*^, unit A'^- vector spins a defined at 
the sites of the lattice, and the Hamiltonian 

rL= -NY,J{x-yyx-(ry-NhY,<yl- (3) 

x,y X 
The partition function is simply 

Z= jX{[da^5{cTl-l)]e-^'^, (4) 

X 

with (3 = 1 /T. We will be interested in studying the finite-size behavior of the theory. 
For this purpose, we consider a finite box Ay of finite extent M in the first q directions 
(called the "parallel" directions and denoted by the subscript ||) and L in the remaining 
p directions (called "transverse" and denoted by _L), with d = q + p, and therefore of 
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volume V = M'^L'P. In order to be able to consider long-range interactions, we define 
a finite-size coupling Jm,l{x) as 

Jm,l{^)= Y1 Y1 J{x + n\iM + n±L), (5) 

and the finite-size Hamiltonian 

n=-N Jm,l{x - y)cT^ ■ (Ty - Nh Y al, (6) 

with periodic boundary conditions. Note that definition implies for the Fourier 
transforms (with p G Ay, see Eq. (fTTH) below) 

MM = E ^''""MLix) = Yl e'P-^Ji^) = JiP)- (7) 

We consider anisotropic long-range interactions and thus we assume that asymptoti- 
cally J{q) has the form 

J(g) ~J(0) + a^ I gj_|2^ + aj| I for|g|^0, (8) 

with < p,a < 1 and a_L < 0, ay < in order to have a ferromagnetic system. For 
simplicity, we will assume the two metric factors to be equal, and, by redefining the 
inverse temperature /3, we can set a± = ay = —1/2. 

In the large- limit, assuming periodic boundary conditions and h = 0, the theory 
is solved in terms of the gap equations 

where K{q) = -2{J{q) - J{0)) and Ay is the lattice 

A*y = (2^M-i ZIj, 2ttL-^ Zl, ) . (10) 

In infinite volume, the same equations holds, with the simple substitution of the sum- 
mation with the normalized integral over the first Brillouin zone [— 7r,7r]'^. 

The meaning of the parameters Ay and ay is clarified by considering the magneti- 
zation and the two-point function. If is the mean value for a system of volume V, 
we define 

Mv={a^)v, Gv{x) = {ao ■cT^)v. (11) 

Then 

M.=a., GH<,) = j^£^, (12) 

where G{q) is the Fourier transform of G(a;). As we show in App. ^ in the limit 
V ^ oo and in the scaling limit Ay — > 0, the correlation function has the form 

Gooix) = ef "''^Goo(a;x/e±,oo, a;||/e||,oo), (13) 

where 
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Therefore, even if the model has long-range correlations in the high-T phase, it is 
sensible to look at ^_l,oo and oo as appropriate typical length scales of the system. 
We will refer to them respectively as transverse and longitudinal correlation lengths. 
They are related by 

eii.oo = (15) 

The critical point is characterized by a vanishing mass gap, i.e. 

^±^00 = = 0' and 

by a vanishing magnetization, fioo = 0, so that the critical temperature is given by (see 
Eq. 0) 

^"^^Z-.,.]^!^^' ^^^^ 

which is finite whenever the effective dimensionality D = p/p + q/cr is greater than 

2 and infinite for D < 2 (the system undergoes a zero-temperature phase transition). 
Moreover, given that we will be interested in finite- volume properties, we have My = 
for all values of /?, so that we can set ay = in the gap equation. 

3 Canonical finite-size scaling 

We do not address here the problem of the definition of a finite-volume correlation 
length and we simply use 

^±,v = i\\y = (17) 

which are the finite- volume analogues of the infinite- volume correlation lengths defined 
above. 

In order to perform the FSS limit we must identify the correct aspect ratio. From 
Eq. (jH)) we see that ~ i.e. L^^ ~ M^" . Therefore, the FSS limit should be 
taken keeping S = M/L^'^^, A = p/a — 1, fixed. In App.EIwe show that, in the limit 
L,M,S^ — > 00 with S and ^±y/L (or, equivalently, y/M) fixed, below the upper 
critical dimension, i.e. for < 4, the gap equation takes a scaling form. Above the 
lower critical dimension, i.e. for D > 2, we obtain 

(4^)p{2-D)/2(^_^^)^p(D-2)^2-D ^ + S), (18) 

where z = {ATT)~p/^{L/^xy)P, (i^ is given in Eq. (|16j) . 

A 1 i np/2p)T{q/2a) f p q\ 

^P,g,P,.- (4^)(p+,)/2p^ r(p/2)r(g/2) V 2p 2a)' ^ 

and Ip,q^p^(j{z-, S) is defined in Add. 1X1 cf. Eq. 

For D < 2 there is no finite-temperature phase transition and the gap equation 
assumes the form 

(4^)p(2-D)/2^^p(D-2)^2-D ^ + j^^^^^^^^^^ s). (20) 

Considering the case D > 2, Eq. (|18|) shows that 

%^ = hi[iP-Pc)LP^''''\S] = 42[(/3-/3c)M'^(^-2),5], (21) 
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that allows us to define a "parallel" exponent z/|| and a "perpendicular" exponent i/± 
as 

Analogously, for the finite- volume susceptibility xv = PGv{0) we obtain 

= A-1 = - fi,)LP^^-''\S] = m2-/^,2[(/3 - /3c)M'^(^-2), 5]. (23) 

The corresponding critical exponent 7 associated with the behavior in infinite volume, 
Xoo ~ iPc — , is correctly identified as 7 = 2piy± = 2avi^\ = 2{D — 2)^^. 

4 Noncanonical finite-size scaling 
4.1 General considerations 

In Sec. Olwe discussed the canonical FSS, obtained by keeping S constant. However, in 
many anisotropic systems the correct exponent A is not known and thus an important 
question is what happens if we consider the FSS limit keeping fixed the ratio 

with 6^0. If one has separately defined parallel and transverse correlation lengths, 
one can immediately identify the correct anisotropy exponent A. It corresponds to the 
value for which £,±y /L and y/M both remain finite in the FSS limit. This uniquely 
defines the correct exponent. This is not surprising, since the correlation lengths satisfy 
Eq. ()15() that essentially defines the correct aspect ratio. On the other hand, one may 
not have access to the correlation lengths but only to some zero-momentum correlation 
function, for instance to the susceptibility. Thus, one may ask what kind of scaling 
behavior, if any, is observed if Ss is kept fixed. In other words, we may ask if we can 
find exponents 7, and t'n so that asymptotically, for L and M going to infinity and 
f3 ^ Pc, we observe a scaling behavior analogous to that defined in Eq. (|23|) . i.e. 

X = L'y/'^f^,i[iP-Pc)L'/'\Ss] (25) 
= M^/^II/^,2[(/?-/?c)mV'>II,55]. (26) 

Of course, the second important question is the relation of these exponents with the 
correct ones defined above. Again, we will restrict our attention to the case 2 < D < 4, 
i.e. below the upper critical dimension and above the lower critical one. 

In the following, we only consider the case (5 > 0, given that the case —1 — A < 6 < 
is obtained by performing the substitutions 

p^q, L^M, p^a, 6 ^ -da^ p{p + aS), Ss ^ SJ^^^^^^'^^K (27) 

If 5 is positive, S goes to infinity as L, M — > 00. We must then study the gap equation 
in that limit, requiring at the same time Xy — > in an arbitrary way. The detailed 
calculation is presented in App.El where we show that there are two nontrivial cases: 

(a) One can take the limit S ^ 00 at 

XyL^^ fixed. In terms of the correlation 
lengths this corresponds to the FSS limit at iA_y/L fixed. At the same time 
^\\y/M ~ {^±y / Ly/'^ goes to zero, i.e. M increases faster than parallel 
correlations. Therefore, the final "effective" geometry is 00'^ x L^, so that we 
obtain the FSS behavior of a "layer." 
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(b) Alternatively, one can take the limit S — > oo, with z ^ 0, keeping z^S"^" ~ XyM'^'^ 
fixed. This corresponds to keeping fixed the ratio y/M while at the same 
time i±y/L diverges. The expected "effective" geometry should be that of a 
q'-dimensional hypercube with linear size M (Af x 0^). 

The first limit always exists and it is easy to predict the result. It correspond to the 
standard FSS limit for a system oo'' x U' and thus the scaling in terms of L is canonical. 
Therefore, it corresponds to taking (5 ^ Pc and L ^ oo keeping fixed (/? — Pc)L^^^~'^\ 
which is indeed the correct combination. Then, we obtain 

X = L2p/^,i[(/3-/3c)i>^(^-'\oo] (28) 

and, using Eq. (j^lj). 



p(D-2) 

(/3-/3c)(^J 



(29) 



Thus i'± = i'± and 7 = 7 (note that both equations give the same result for 7). On 
the other hand, the parallel exponent z^n is given by 

l + A + 5 

a result which is a trivial consequence of the fact that if Ss is fixed, then i'\\/i'± = 
1 + A + 6. Therefore, even if 5 ^ 0, one still observes scaling. Transverse exponents are 
correct, while longitudinal ones are only effective. Let us finally note that Eqs. (|28() 
and H29|) give both the correct infinite-volume behavior for x- Indeed, the existence of 
a finite infinite-volume limit for j3 < j3c implies fy^^i{x) ~ (— 2;)"^/*^^"^^ ~ {—x)~'^ for 
X — > —00 and therefore Eqs. (PH|) and 1)29(1 are consistent with Xoo(/?) ~ (/?c — ■ In 
other words, one must have 7 = 7 in order to have the correct infinite-volume limit. 

Limit (b) is much less conventional and we will show that such a limit exists (in the 
sense that x has a scaling behavior) only if q/a < 2, i.e. if the g-dimensional theory is 
below the lower critical dimension: in this case however, the exponents obtained from 
the data collapse do not have anything to do with the correct ones, while the scaling 
functions depend only on the ^-dimensional theory, i.e. theories with different p and/or 
p — they are thus physically inequivalent — have the same scaling functions. Thus, in 
this case from the observation of a good scaling behavior one can draw the incorrect 
conclusion that physically inequivalent theories have the same critical behavior. 

At this point the reader may be puzzled by the fact that for the same Ss we observe 
two different types of scaling. Mathematically this is related to the fact that in the 
FSS limit we must follow a path in the /?, L plane with (3 ^ (3c and L — > cxd. The two 
limits correspond to two different families of paths. Unexpectedly scaling is observed 
on both of them. 



4.2 Limit (b) 

We wish now to consider the limit ^ 00, z ^ at fixed zS" . Such a limit has 
been considered in App.O where it is shown that the leading contribution to the gap 
equation has the form 
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where r = {A-n)p/"-^ S"^ , uj = zr"!'^ = (47r)-'^/2Aj/^M'^, Ag^^ = ^o,g,p,a (it is easy to 
see that it is p- independent), and Ig °{uj) is the analogue of Ip^q^p^a{z, S) appearing in 
Eqs. H18() and 1)20^ for an isotropic (^-dimensional system with long-range exponent a. 
The constant K, cf. Eq. (|8.S)) . gives a subleading correction for q/a < 2. We will now 
distinguish two cases, depending on the value of q/a. 

4.2.1 Case < g/a < 2 

In this case we can neglect K and set r = rsL^^ , where rs = {4:Tr)P^'^~^ Sg (by definition 
ro = r) is kept fixed in the limit. Then, for /3 > /3c, we can rewrite 

(4^)p(2-D)/2+p/2^./2-a(^ - p^)LP(D-2)-2Sia-c/2)^2~,/. ^ j^so^^^ _ ^^^^^ (32) 

This equation has the same functional form as the gap equation of an isotropic q- 
dimensional system with effective dimension des = ^/c < 2, cf. Eq. (|20)) . Therefore, if 
for a g-dimensional isotropic system of extent M and long-range exponent a we have, 

then in the limit we are considering we obtain 

X = M^^f^°[S^''''^^+''+'^M^''\\ (/? - (34) 

where 

^ (1 + A)z.|7i + 6a{d,s - 2) 
II l + A + (5 ^ ^ 

We note that the second term in the numerator of Eq. (|35j) is equal to ^^^^1^^, where 
^||,iso = [o"((iefr — 2)]""*^ is the correct exponent for the dgfr-dimensional isotropic system, 
so that Eq. (|H5|) interpolates between and ^^l^o- 

There are two limitations to the validity of Eq. (|H4j) . First, Eq. (|34j) holds only 
in the limit in which the argument is constant as M ^ oo. Since we are considering 
/? — > /?c, this requires t'n > and in turn 

(D-2)(l + A) 

5 < 5.nax = (2 _ rf^^) • (36) 

If 5 is larger than Jmax no scaling is observed by taking (3 ^ (3c- One observes a truly 
Q'-dimensional FSS behavior: a scaling behavior is obtained only by taking (3 — > oo. 

Additionally, it should be observed that f^°{x) is defined only for x > 0, and indeed 
a solution to Eq. H32() is found only when the left-hand side is positive. In Appendix 
Owe show that the correct extension for x < is f^°{x) = in the sense that, for 
P < /3c, we have x/M^"^ ^ 0. 

We also obtain 

7 = 2crPj|, 

The effective exponents are thus unrelated to the correct ones and vary continuously 
with 6, interpolating between the correct ones {S = 0) and those of a deff-dimensional 
system {5 = (5max)- 
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At this point the reader may be puzzled by the fact that "i, since we already 
claimed in Sec. 14.11 that this is a necessary condition to ensure the correct infinite- 
volume limit. We now show that there is no contradiction. The effective exponent 7 
differs from 7 because in limit (b) the scaling function f^{x) vanishes identically in 
the whole high-temperature phase x < 0, i.e. in the only phase in which x is finite 
(Xoo(/3) = 00 in the whole low-temperature phase). In other words, scaling (b) does 
not give any information on the infinite- volume behavior, and thus there is no surprise 
for 7 7^ 7. 

4.2.2 Case q/a > 2 

In this case K in Eq. cannot be neglected, and therefore no scaling solution can be 
found. Scaling can only be observed along a very specific trajectory in the P, L plane. 
Indeed, if we consider the line 

=/3, + aL-i/^^, (38) 

where 

a = ^4^yp{2-D)/2-d/2j^^ (39) 
we still observe scaling in the sense that asymptotically 

X = - /3,(L))], (40) 

where uj^^ = uj^^ + 6a{dc,s — 2) = uj^^ + ^^\\iso fx° scaling function of a 

deff-dimensional isotropic system (with defr = q/cr > 2). However, this is not standard 
scaling since a must be tuned to a proper value: it is not enough to take /3 — > /3c, but 
it is also necessary to take this limit along the line ()38|) . 



5 A numerical example 

We wish now to present a numerical example, in order to clarify the issues discussed 
in the previous sections. We consider the standard short-range model with nearest- 
neighbor couplings, so that p = a = 1 and A = 0. We consider a three-dimensional 
cubic lattice of dimensions M x L^, i.e. p = 2 and q = 1. For this system \A2\ 

a/3 — 1 

13c = ^r(l/24)2r(ll/24)2 0.25273100985866300303 . (41) 

1927r'^ 

According to the results of Sec. 14. 21 we should be able to observe a noncanonical scaling 
for < 5 < 1. We thus fix (5 = 1/2 and compute x for several values of f5, L, and 
M corresponding to S'1/2 = 1 and S1/2 = 10. We use 5 < L < 245 and 0.235 < /3 < 
0.279. The previous results predict two possible scalings. Scaling (a) corresponds to 
consider x/M^^^ vs (/? — f3c)M'^^^, while scaling (b) corresponds to consider x/-^^ vs. 
(/3-/3c)mV3. 

In Fig. H we report the results for limit (a). For S'1/2 = 1 there are somewhat large 
corrections to scaling and the convergence is slow in the low-temperature region. For 
'S'1/2 = 10 tfis agreement is very good: no corrections can be seen on this scale. This 
different behavior is due to the different lattice sizes that are used: the data have the 
same values of L and therefore lattices with S1/2 = 10 are ten times larger than those 
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Figure 1: Noncanonical FSS scaling in 3 dimensions: limit (a). The solid line is the theoret- 
ical prediction. Here S = 1/2. 

with Si/2 = 1- In Fig. [2 we show the same data in a log-log plot. Here we see better 
the results corresponding to the high-temperature phase, although corrections in the 
low-temperature phase are somewhat less visible. 

In Fig. we report the results for limit (b). The data show a reasonably good 
scaling behavior, although the convergence to the limiting curve, corresponding to a 
one-dimensional system, is quite slow. For limit (b) we have also investigated the 
behavior of the leading correction term. In App. IC.3I we compute the first scaling 
correction to Eq. (|3l|) . It can be effectively taken into account by a size-dependent 
shift of the critical temperature, i.e. by writing 

XM~^ = f^°[SJ^^\p - fic)M^'^ - (42) 

where K = K{4:7r)~^S^ K is defined in Eq. (|83() . Numerically, in this particular 
case, K -3.9002649200019558828. In Fig. H we plot the difference between the 
left-hand and the right-hand side of Eq. H42|) versus the argument of the function 
/3f°. It is clear that data points in Fig. |1] are converging towards zero, confirming 
that asymptotically f^°{x) is indeed the correct scaling function. A numerical analysis 
shows that corrections to Eq. ()42() vanish as M~^. 

From a practical point of view it is important to be able to distinguish scaling (a) 
from scaling (b). This is not very difficult since in scaling (b) fx{x) = in the region 
x < 0. This observation provides a simple criterion: if the high-temperature data scale 
onto a nontrivial curve we are clearly observing limit (a): for instance the scaling curve 
appearing in Fig. |2l can only correspond to scaling (a) . 

6 Conclusions 

In this paper we have considered an anisotropic O(cxo) model, focusing on the shape 
dependence of the FSS limit and in particular on the role of the aspect ratio. Consid- 
ering for instance the susceptibility, we have shown that, even if the FSS limit is taken 
by keeping constant an incorrect aspect ratio M / L^~^^~^^ , one still observes scaling, 
i.e. the data satisfy the scaling forms (|25|) and (|26j) . although the effective exponents 
differ from the correct ones. Therefore, the numerical observation of data collapse does 
not give direct information on the correct values of the anisotropy exponent A and 
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Figure 2: Noncanonical FSS scaling in 3 dimensions: limit (a). The same data of Fig. in 
a log-log plot. 




Figure 3: Noncanonical FSS scaling in 3 dimensions: limit (b). The solid line is the theoret- 
ical prediction. Here 6 = 1/2. 
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Figure 4: Corrections to the FSS behavior: hmit (b). Here S = 1/2. 

of the critical exponents. In the case of hmit (a), if we consider the scahng at hxed 
M/L", we obtain for a > 1 + A 

z>ll = ah'± 

7 = 7, (43) 

while for a < 1 + A 

i'± = y\\la 
= 

7 = 7. (44) 

These expressions have been explicitly verified in the specific model we consider, but 
we expect them to be valid for generic anisotropic systems. In the case of limit (b) 
we are not much interested in the specific form of the effective exponents. The only 
important feature that we must notice is that all exponents vary continuously with a 
and, if a is sufficiently different from 1 + A, scaling is observed only for /3 — > oo. 

From a practical point of view, it is very important to establish a strategy that 
can be used to correctly identify the exponents in anisotropic systems. One possibility 
consists in performing a numerical study keeping fixed MjU^ for different values of a. 
If one observes that the data obtained with two different values of a collapse onto a 
single curve by using the same exponents and the same value of /3c, then the exponents 
which are determined in this way are correct. Indeed, this implies that one is observing 
limit (a) — in limit (b) all exponents vary continuously — and therefore Eqs. (|43|) (|44|) 
apply. Note that as an additional check one can see whether fy^{x) is nontrivial for 
X < 0: if this the case, one is observing limit (a). To be more specific, suppose that 
the scaling behavior ()25|) is observed for ai and 02, 7 and i'± being the same in both 
cases; then 7 and i'± can be identified with 7 and i'± and 1 + A < min (ai, 02)- From 
the data at ai and 02 one cannot determine i^y and A. For this purpose, one must 
perform simulations at a value as such that u± 7^ i.e. one must find a value such 
that as < 1 + A. Then, one should use Eq. setting 7 = 7 (this is meant to avoid 
limit (b)) and using as a free parameter. Then z^y = and 1 + A = i/||/z/_|_. As a 
check, one can analyze the data at ai and 02 and look for the scaling behavior (PH]). 
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fixing 7 to avoid limit (b) . The effective exponents should satisfy = az/y /(I + A) . 
Thus, in principle, by studying the behavior for three different values of a one can 
determine all exponents. However, in practice this may be difficult, because it requires 
the ability to distinguish corrections which vanish as M, L — > oo from corrections that 
persist in the limit M, L — > oo and that require therefore the use of different values of 
the effective exponents. 

From a practical point of view it is important to note that all the above-discussed 
ambiguities are not present if the correlation lengths are measured, since their scaling 
behavior is fixed, ~ L, ~ M. Deviations from these laws give immediately 
the exponent A. Therefore, finite-size numerical studies of anisotropic systems should 
always determine finite-volume correlation lengths, as done for instance in Ref. [33] for 
a driven lattice gas. 

A somewhat unexpected result of our analysis is that in some cases it is possible 
to observe two different types of effective scaling: one corresponding to the layer /strip 
geometry oo"^ x L^, which is the scaling one would naturally expect, and the other one 
corresponding to a lower-dimensional system M'^ x 0^. The possibility of this type 
of scaling is not obvious since it requires that one correlation length increases much 
faster than the shortest of the two lattice sizes. Even less clear is why this limit exists 
only if q is below the lower critical dimension of the model (for larger values of q such 
a scaling can still be observed by performing an additional tuning). We expect that 
scaling (b) is strictly related to the fact that x ^ as defined here diverge in the 
whole low-temperature region, but it is difficult to transform this conjecture into a 
quantitative argument. In Ising systems it is simple to avoid this feature. For instance, 
one could define a finite- volume susceptibility as 



An open question is whether such a quantity may show a noncanonical scaling (b). 
Our guess is that only scaling (a) is possible in this case. 

A Gap equation: Finite-size scaling limit 

In this Appendix we wish to determine the asymptotic form of the gap equation for 
Xv 0, M,L oo, without making any hypothesis on the way in which the limit 
is taken. We will assume to be below the upper critical dimension, i.e. L) < 4 where 
D = p/ p + q/a. We will consider separately the cases D < 2 and D > 2. 

A.l Gap equation above the lower critical dimension: 2 < 




(45) 



where 




(46) 



L> < 4 



We start by rewriting (see Ref. j^H]) the gap equation as 



/3 - /3e = -Ay 




d'^p 1 



(47) 



{27rrK{p){K{p) + Xv) 
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where 

1 ^ 1 r d'^p 1 

= V ^j^^ K{p) + Xy ~ y[„,,,]. {2^YK{p) + \v- ^ 

Now, for D < 4 and Ay ^ we can expand K{p) in the integral appearing in Eq. (|47|) 
and extend the integration over all W^, obtaining 

/3 - /3e « -Xy'^'^Bi + S, (49) 

where 

d'^p 1 



with Kc{p) = + and ^p,g,p,o- defined in Eq. (|T9]) . Using the Poisson 

summation formula for a periodic function /(x) = f[x + I), 



l—l oo „; 



(51) 



n=0 k=~oo ' 

we can rewrite 

^ (2vr)^ K(p) + Av ^I>^(2vr)'^ Kc(p) + Ay ' ^^^^ 

where in the last step we have assumed L,M ^ oo. Convergence is guaranteed by the 
oscillating phase factor. 

In order to compute the asymptotic behavior of S for Ay — > 0, M, L oo, we 
introduce the Laplace transform /^(p) of exp(— -u"), i.e. we define 

'dpe-'^PUp) = e-^", (53) 

for < a < 1. Of course, for a = 1, fi{p) = 5{p — 1). It is easy to show that for 
< a < 1 

/a(p) ~ forp^+oo, (54) 

and 

fa{p) ~ p-(2-«)/[2(i-")] exp [-(1 - a)(p/a)-"/(^-")l , for p ^ 0+. (55) 
Then, we obtain 

{'^^Y Jo 

fll 7o JK'' (271")'' Jo 



M 



= (4^)p+(d-pD)/2 d-qdr fp{r])f^{T)Hp^q{r],T,z,r) = Xy' Ip^g^p^„{z, S), (56) 
where 5 = M/Lp'", r = {Att)p/''^^ S"^ , = {4^)-pl'^p\v, Vt = t^^'^V, n = t^^'^r, 

Hp,,{v, T, z,r)^ dt ry-^/Vr^/'e-^'* [B^ {^i') B'^ (rrr^) - l] , (57) 
Jo 
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and 

Note the well-known property 

B{s) = s'^^^B{l/s), (59) 

that implies B{s) ~ s~^/^ for s — > 0. 

Using Eqs. (|i9|) and (|56|) we obtain Eq. (fT8|) . Note that for an isotropic geometry 
with g = 0, Ipfl,pfi{z, S) does not depend on S and thus we will simply write Ip°p{z) = 
Ipfl,p,oiz, S). For a geometry with M = oo — we call it layer geometry — Eq. (fT8|) still 
holds with S = oo. In Eq. (|57() it corresponds to setting r = cxd, i.e. replacing B'^irr^^) 
with 1. 

In the low-temperature phase (/3 > /3c fixed) we have Ay for all (3 when 
V — > oo. Thus, we can use the above-reported expressions to determine the infinite- 
volume behavior of ^j__y. From Eq. (|49|) we obtain S = (3 — 13c > in the infinite- volume 
limit. Because of the prefactor L^'^~'^)f appearing in the last term of Eq. ()56() . since 
D > 2, this is possible only if the integral, i.e. Hp^g{r],T, z,r) diverges (the r/ and 
T integrations are completely harmless), which only happens for z — > 0. In order 
to compute the asymptotic behavior of Hp^q{r], r, z, r) for z ^ 0, we use the duality 
property ((S^ to rewrite 

HpMr,z,r) = [\tr^;'/\-'^/'e-''[B^^^')B'^{rri')-l] 
Jo 

/oo 
dtr-^/^e-'"' [BP MB'i (rt/r) - 1] 

/"OO , , 

+ / dt e--'* (r-«/2 - . (60) 



The first two integrals are always finite while the third one gives Hp^q{r],T, z,r) 
T~il'^ j z^ . This implies that in the infinite- volume limit 



A. 2 The case D <2 



For D < 2 little changes. Now criticality is observed only for /3 — > oo. We rewrite the 
gap equation as 

f d'^p 1 
^ = W (2vr)^K(p)+Ay+^- ^''^ 

For D < 2 and Ay — > 0, we can expand K{p) in the first term and extend the 
integration over all W^, while S can be treated as before. We obtain 

/? = \^'^'^B2 + A{?/'"^/p,,,p,.(z, 5), (63) 

where 



B2 ^ I , . = -^P,9,P,<x, (64) 



d'^p 1 
(27r)rf Kc{p) + 1 

and Ap q^p^fj is defined in Eq. (|19|). Eq. (|20|) follows immediately. 
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B FSS of the correlation function 

Now we consider the (real-space) correlation function Gv{x): 

Gv-(-) = f E j^^- CSS) 

and compute its asymptotic behavior for L,M,\x\ — > oo and Ay 0. Using the 
Poisson summation formula (|5H) we obtain 

Gv{x)= r dte-^^'Y ! .^^ip\l{n^^M+x^O+ipAn^L+x^)~tKip) ^ ^gg^ 

^0 „ Jl-7T,7r]'i (2vr)'' 

Since we are interested in the limit \x\,L,M oo we can replace K{p) with its small-p 
expansion and extend the integration to the whole space, obtaining 

By rescaling t, we obtain finally 

^ g-|(xx/L)+nxlV(4r?t)-.S2|(^||/J\^)+n||lV(4rt) 



Lp'^^"D)F{x±/L, xu /M, S, XlL^p) 



(68) 



C Gap equation for 5* ^ oo 
C.l General results 

Now we will determine the limiting form of the gap equation (|47|) for L ^ cxd, A/ ^ oo, 
S ^ CO and arbitrary values of z = ^^"^ {L / ^±yY , including z = 0. Indeed, we 

shall mainly be interested in the limit z ^ together with ^ oo. 

For this purpose we need to compute the asymptotic behavior of Hp^q{r], r, z, r) (see 
Eq. (|57j) ) for r — > oo. We start from Eq. (|6())) . For what concerns the first term in 
this equation, we can take the limit naively, i.e. we can replace B{rTf^) with 1. The 
second term can be rewritten in the form 

/oo rco 
dte~'''[BP{r]t)B'^{Tt/r) - 1] = r''^/^ J dt e~ [B^ {r]t) - l]B''{Tt/r) 

dt e-''^B'i{Tt/r) - r-i/^ J dt e"^"*. (69) 

In the first term we can simply replace B{Tt/r) with its asymptotic behavior for r^/r 
small, i.e. {Tt/r)~^/^ , while some additional manipulations are need for the second one. 
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Rescaling t = sr'^, we rewrite it as 

/oo rco 

Jl ^ 

In the second term, one can easily convince himself that it is safe to replace the lower 
integration limit r~'^ with 0. Collecting everything together, we obtain 

Jo 

/oo 
dte-'\-'i'^'' [BP{r]t) - 1] 

dte-' [B'^iTt) - 1] 
+r--''/2 I' dte-^'^^'r;''''' [S" (r,"') - l] 

/>oo 



Then, for D > 2, we can write the gap equation as 
where 

Mx) = ^ + (73) 

/oo 
d«-^/2e--* + g(2),,,,(x), 



(74) 



and 



/•oo 

gW(x) ^ / dT/,(T)GW(T,x), 
JO 

GW(t,x) ^ l^^dte-^''[D"{n)-l]+j\te-^'T-'^/'[B'^{Tr')-l] 

= J dte-^V-«/'[E«(T-i)-l]+r-«/2y die-^*r«/2^-^, 

/•oo 

OSL,i^) ^ C,,J dnfp{^)Gfl,^^{^,x), (75) 

^0 

Gg,,,,(r/, x) = dte-^H-'^l^'^ [BP {m) - M + dte-^H-'^/'V^' [B^ {rj^') - l] • 
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Note that Q^l{x), Qj^l^a,q{x), and J2{x) are finite for x — > and the relation 

/oo 

/•oo 

-T-«/2 / dte-'\~'^/^''. (76) 



We wiU be interested in computing the behavior of J7i(x) for x ^ oo. If < a < 1, we 
rewrite 

Ji{x) = / drUr) / dte-^'r;'''^ {r^^) - l] + C,,, / dt e-'^H-'^/^^ (77) 
Jo Jo Jl 

By performing an integration by parts it is easy to see that the last term vanishes as 
e~^/x. The asymptotic behavior of the first term depends on a. For a = 1 we can 
rewrite it 

/ dte-'''H~'i/^[B'i[t-^) -I] ^ 2q t-'?/2g-at-^t-i 
Jo Jo 

/•oo 

Jo 

/7r\3/4-g/4 

« 2q(-) (78) 



For < cj < 1, more work is needed. First, we replace r with s = rt^^'^ = Tt obtaining 



fOO fOO 

/ dt t-i/-e-^* / ds s-'^'^Usr^'n [B'i {s-^) - 1] . 
Jo Jo 



(79) 



In the limit of large x, only the small-t region contributes to the integral (because of the 
exponential factor e~^*) and the function fa{st~^f'^) can be replaced with its small-t 
behavior (this is safe because the integral over s is convergent near s = 0), i.e. we can 
replace fa{x) with x^^^'^ . Thus, the integral (|79j) becomes 

poo poo 

/ dt t e-"* / ds s-^-^-"/^ [Bi (s-i) - 1] 
Jo Jo 

= dss~'~--^/^B^{s-')-l], (80) 

and it is easy to check that for a > 0, q > the second integral is convergent. Therefore, 
as X — > oo, we have that J7i(x) ~ x^^ for any < a < 1 and J^i{x) ~ x^ exp(— C-^/x) 
for (T = 1. 

From Eq. ()72() we see that there are two interesting scaling limits: (a) r — oo at 
fixed z; (b) r — > oo at fixed z'^r'^. As we will show below these two limits correspond to 
the two different cases we discussed in Sec. HI They are descussed in detailed below. 



C.2 Limit (a) 

In this case we should consider the limit r ^ oo at fixed z and we should there- 
fore reobtain Eq. (|18)) with 5 = 00. This follows immediately from the fact that 
r''-i/'^Ji{z^r'') 0, cf. Sec. EH and from 

/OO roo 
(ite-^V«/2-= d7]dTfp{7])UT)Hp,,{v,T,Z,00), (81) 

which is a direct consequence of Eq. (|76|) . 
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C.3 Limit (b) 

Now we discuss the limit ^ and r — > oo at constant z^r"^ . Note that the last term 
in Eq. (|72|) behaves differently, depending on whether q/2(j is larger or smaller than 1. 
For q/a > 2, integrating by parts, we obtain 



2a __,2 
-e ^ 



q-2a 
2a 



2a 



(82) 



Then, note that for D < 4 we can extend the remaining integration down to zero. 
Moreover, we can approximate with 1 since is subleading with respect to 

^o--g/2. ii^jggjj ^2 ^ ^-cr ^ ^a-<3/2 ^ ^g/2-2cr q < Aa foT D < 4. Therefore, if we 

define 



K = J2{0) + 
Ji{x) = Ji{x)- 



2a 



q-2a 
2a 



-a 



(83) 



2a 



Ccrj 



/•CO 

= / drdt^(T)e~^Vr^/'[i3'^(T,-^) - 1] + a,,r(l - q/2a)x'-''-^'^^/^%{8A) 
Jo 

we obtain 

where we have also discarded the term proportional to z^^^ ~ {^^^(jl2\p-2^(j~ql2^-cTpl2p ^ 
which is subleading with respect to r'^~il'^ . 
For q/a < 2, we write 



-z^t^-q/2a 



a-q/2 / ^-zh^^t^-q/2a 



a^g/2 / ^-z\'^t^~g/2a 



Jo 



It is then easy to see that we reobtain Eq. (|85j) . In this case K is subleading and can 
be neglected in the scaling limit. 

In this derivation we have implicitly assumed that it is possible to take the limit 
z ^ 0, r — > oo, with z^r" fixed. However, we now show that this is not the case for 
/? < fic- Indeed, the function J\[z) is positive and thus Eq. can only be valid if 
/3 > /9c- If this condition is not satisfied, the only possibility is that (/? — Pc)L^^^~'^^ 
vanishes in the scaling limit as well as the right-hand side, which in turn implies 
z^r"" oo. This means AyM^"" oo, or x/M'^'^ 0. 

Finally, we rewrite the gap equation in terms of the constant Aq^^y = ^o,q,p,o- (it is 
easy to check that it is p-independent) and of the integral Pq^'^{z) that appear in the 
gap equation of an isotropic g-dimensional system. It is indeed easy to show that 



Jl{z^) = {A7Ty/^Z-'+'^/^ [/^°(Z) - Aq,,] . 



(87) 
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